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We propose a universal quantum computing scheme in which the orthogonal qubit states |0) and 
|1) are identical in their single-particle spin and charge properties. Each qubit is contained in a 
single quantum dot and gate operations are induced all-electrically by changes in the confinement 
potential. Within the computational space, these qubits are robust against environmental influences 
that couple to the system through single-particle channels. Due to the identical spin and charge 
properties of the |0), |1) states, the lowest-order relaxation and decoherence rates 1/Xi and I/T2, 
within the Born-Markov approximation, both vanish for a large class of environmental couplings. 
We give explicit pulse sequences for a universal set of gates (phase, n/8, Hadamard, CNOt) and 
discuss state preparation, manipulation, and detection. 

PACS numbers: 03.67.Lx, 73. 21. La, 85.35.Be 



X 



I. INTRODUCTION 

Proposals for quantum computing architectures based 
on semiconductor devicesi^^aj££ are attractive for 
their scalability; once the few-qubit problem is solved, 
massive scalability is not expected to pose insurmount- 
able barriers either in resource requirements or fabrica- 
tion precision. This is primarily due to the sustained 
and continued improvements in epitaxy and lithography 
over the past few decades, and the ability with which 
new techniques, often developed in industry, are trans- 
ferred to basic research laboratories. On the other hand, 
semiconductor environments are hardly systems of pris- 
tine quality and isolation, and there are severe trade-offs 
between long coherence times and short access times. 

Pure spin qubits, for example, couple relatively weakly 
to their environment & Their dipole tails are often neg- 
ligibly weak and spin exchange effects, while potentially 
strong, are short range. But precisely because of this 
weak environmental coupling, spin qubits may be po- 
tentially difficult to control and manipulate. For single- 
particle qubits, local Zeeman tuning is required to rotate 
bits. The opposite scenario is often true for charge qubits. 
Here, control may be attained very quickly with metallic 
gates or optics. 6,7 However relaxation and decoherence 
times can be very fast, requiring even faster switching 
times. 

This begs the question of whether there exist hybrid 
qubits which accentuate the positives and mitigate the 
negatives. We show below that this does indeed seem 
the case if a single qubit is judiciously defined as a corre- 
lated few-body system whose charge and spin degrees of 
freedom are entangled. These strong correlations should 
additionally be effective at suppressing relaxation and 
decoherence through single-particle channels. Indeed our 
two orthogonal qubit states are identical in their single- 



particle spin and charge degrees of freedom; differences 
only show up in their two-body correlation functions. 

Sources of decoherence and dissipation can also be 
broadly classified as spin based or charge based. Both de- 
stroy the unitary dynamics of the system either by taking 
it outside the computational subspace, or by remaining 
within the computational subspace, but causing either 
uncontrolled qubit flips, or pure dephasing without dis- 
sipation of energy. This will generally occur whenever 
an environmental influence couples differently to each 
qubit state. For example, if the two qubit states dif- 
fer in their spin, then random magnetic fields are an is- 
sue. For single-particle qubits, this will always be the 
case, and likewise for two-particle qubits; it is not pos- 
sible to define two orthogonal one- or two-particle states 
with identical spin and charge densities. A three particle 
system, however, can be constructed in which both the 
charge density and the spin of the two orthogonal |0) and 
|1) states are identical. 

We show that the qubits we define below admit a uni- 
versal set of one and two-qubit gates, and we give ex- 
plicit gate pulse sequences which implement this univer- 
sal set. We also discuss issues of decoherence and relax- 
ation among the qubits and show that, for a broad class 
of environments, including certain spin dependent ones, 
relaxation and dephasing are absent (1/Ti = 1/T V = 0) 
within the lowest-order Born-Markov approximation. We 
expect the residual decoherence rate due to higher-order 
couplings, non-Markovian effects, and other, weakly cou- 
pled, environments to be small. We also discuss exten- 
sions to the model of system-environment coupling, and 
comment on issues of state preparation and detection. 

In the following section, we describe our model 
electronic Hamiltonian consisting of two many-body 
parabolic-elliptic quantum dots, with long-range in- 
tradot Coulomb repulsion. In Sec. IHII we construct our 
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qubits and demonstrate how correlations produce orthog- 
onal |0) and |1) states with identical spin and single- 
particle charge densities. Section |IV] contains explicit 
implementations — in the form of pulse sequences — for a 
universal set of quantum logic gates (Hadamard, 7r/8, 
Phase, and CNOT gates). Section fVl demonstrates that, 
to lowest order, intra-qubit relaxation and dephasing is 
absent for all pure spin and pure charge environments 
which couple to the qubit through single-particle chan- 
nels. Finally, in Sec. IVIl we briefly discuss issues of state 
preparation and detection. 

II. MODEL HAMILTONIAN 



The long-range Coulomb interaction Hq ou i can then 
be written in the usual second-quantized form as 

-ffCoul = ^ ^2 V Vklct a cl a ,Cl a >C ka , (5) 

where all indices are summed over. An explicit and exact 
closed form expression for the matrix element 

f e 2 
J 27rge 

_ (6) 
is derived in Ref . [lfj . 



We consider two coupled dots with Hamiltonian 

H = -ffdotl + -ffdot2 + ffcoupl, 



(1) 



where the first two terms denote the individual quantum 
dots whereas the third denotes interdot coupling. For 
-ffcoupi, we shall take a simple coupling Hamiltonian but, 
within each quantum dot, we shall take full long-range 
repulsive interactions into account (exactly). We first 
focus on a single qubit and subsequently discuss two- 
qubit interactions. 

We encode a single qubit in a single elliptically 
confined^ two-dimensional lateral quantum dot. We 
place three interacting electrons in the dot and con- 
sider the two-dimensional subspace spanned by the S — 
1/2, S z = —1/2 spin sector. Single qubit rotations are 
created by tuning the eccentricity of the elliptic confine- 
ment potential^ whereas two-qubit operations, as we 
show below, are created by controlling the coupling be- 
tween two adjacent quantum dotspii 2 - 

The Hamiltonian of a single dot is given by -f/doti = 
-ffibody + -ffcoui, where 
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We take a magnetic field B = (0, 0, B) perpendicular 
to the plane of the dot. The Hamiltonian @ can be 
exactly diagonalized with canonical Bose operators a\, 
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and their Hermitian conjugates a: 
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III. QUBIT CONSTRUCTION 

For definiteness, we consider three singly-occupied or- 
bitals \nm) = |00), |01), 1 02) corresponding to the three 
lowest-energy orbitals in the lowest Landau level. With 
this orbital occupation, the S — 1/2, S z = —1/2 sub- 
space is two-dimensional and is spanned by our orthogo- 
nal qubit states 

|0) = -^(2|UT)-|iU>-|TU)), (7a) 



|i) = _(|iti)-|TU». 



(7b) 



Each term on the right is a single antisymmetrized state: 
IS0S1S2) = c oos c oisi c 02s 2 l vacuum )> where the operator 
c lims creates an electron in state \nms). The spin con- 
figurations in Eq. |J7J| are (up to an overall exchange of 
spin up and down) those in Ref. [12]; however, the elec- 
tron states differ in their orbital degrees of freedom. In 
particular, the states (jTj) have particles sitting in orthog- 
onal orbitals; this orthogonality is required for the charge 
densities to be identical during gate operations. 

The states in Eq. (|7J) cannot be written as single Slater 
determinants in any single-particle basis; they are corre- 
lated states with entangled spin and charge degrees of 
freedom. These correlations enable the states to be both 
orthogonal to each other and yet exhibit identical sin- 
gle particle properties. Both qubit states in ([7]) have 
spin S — 1/2, S z = —1/2. Furthermore, defining the 
charge density operator as p{r) = ^ 8{r — fj), we find 

(p)\o) = (p)\i) = Ei IV'oiMI 2 , where Vw(r) = (r\nm) is 
a real-space eigenstate of Eq. © . The density is plotted 
in Fig. [T] for two different values of z = lu c /lo x . 

Physical differences in the qubits arise at the two-body 
level. For the two-particle density 
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We build many-body states through antisym- 
metrized products of single-particle states |nm), 
where a\a\\nm) — n\nm) and a\a2\nm) = m\nm). 



we find 



Sp — 2Fqi — Fq2 — F12, 



(9) 
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FIG. 1: (Color online) Identical charge density (p(r)) for both 
qubit states. Both plots have uj y /uo x = 1/2. The left plot is 
at zero magnetic field whereas the right has u) c fuj x = 5. For 
Ux = 1 meV, this corresponds to B z ~ 3 T for GaAs. 
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FIG. 2: (Color online) Two-point density (pt(ri, r<i))\Q) with 
ri =0. The left plot is for Q = and the right for Q = 1. 
Both plots have u> y /cu x — 1/2 and lu c /lo x = 5. 



where (5p = (/?t 



(ft 



|o) 



and 



2^ = Re [^(r 1 )^(r 2 )^(r 2 )V*(r 1 )] . (10) 

The two-point functions are shown in Fig. [5] for r\ = 0. 

Because both the single-particle charge and spin prop- 
erties for both qubit states are identical, intra-qubit de- 
coherence and dissipation should be minimized. We show 
below that, within the lowest-order Born and Markov ap- 
proximations, the T\ and T 2 times are infinite for a very 
large class of environmental models. Before doing so, 
however, we first show that a complete universal set of 
logic gates is achievable in this system. 



IV. UNIVERSAL QUANTUM LOGIC GATES 

In the space defined by the qubit states in Eq. ([7]), 
the electronic Hamiltonian ^ can be written as a 
pseudospin-1/2 particle in a pseudomagnetic field. In 
this particular case, we havei£ 



-ffqubit = b x a x + b z a z + &0<70, 



(11) 



with a x , a z the Pauli spin matrices and &o the identity 
matrix. The pseudomagnetic field components b x , b z , 
and bo are given by 

b x = V3(y 220 -^22i)/2, (12a) 

b z = -Vouo + (Vi 2 21 + V0220) A (12b) 

bo = Voioi + V0202 + V1212, (12c) 



where Vijji are exchange (and Vijij direct) matrix ele- 
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ments, given in Eq. 
m 



with rii = rij = 0, rrii 



and 



— j. Explicit (exact, analytic) expressions of these 
are given in Refs. [Iolll3T ]. 

The main point with regard to qubit rotations is that 
the fields in Eq. (fT2")l have a different functional de- 
pendence on the dimensionless ratios r = u> y /u> x and 
z = lo c /lj x . Thus, adiabatically controlling cither r{t) 
or z(t) can rotate qubits^ These ratios may be changed 
at fixed magnetic field (lo c ) by altering the two confine- 
ment frequencies ui x and oj v independently. 

To perform an arbitrary computation, we require a uni- 
versal set of quantum logic gates which typically consists 
of both single and double qubit operations. We focus first 
on the single-qubit portion of this universal set, followed 
by the two-qubit portion, the CNOT gate. 



A. Single Qubit Gates 

A universal set^ of quantum logic gates is given by the 
CNOT gate, which we discuss below, and the single-qubit 
Hadamard gate H, tt/8 gate T, and phase gate S. These 
are each given by 
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1 1 

1 -1 
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r/4 



s 



(13) 

We expect that all SU(2) operations on qubits encoded 
as in Eq. ([7]) can be achieved asymptotically since two 
non-parallel pseudofields are achievable with two differ- 
ent values of r and 

In order to find explicit time-dependent parameters 
r(t) and z(t) for which the single-dot time evolution 



U = T cxp -i / H dotl (r(t), z{t)) dt 



(14) 



equals the desired single-qubit operation, we adapt the 
minimization method used in Refs. [Hl|]. 
time-ordering operator.) The time interval 



(T is the 
[0,T] is di- 
vided into N discrete pieces during which the functions 
r(t) and z(t) are set constant. We are then left with 
an optimization problem with 3A^ variables U, ri, and 
Zi (i — 1,..,N) where U denotes the length of the ith 
phase, Ylj—i ti — T; the parameters and Zi determine 
the values of r(t) and z(t) in the ith phase. We nu- 
merically minimize the function / = | \U({ti, Ti, Zi}) — 
Ut\\ where U({ti, ri, z{\) is obtained by exponentiation, 
Eq. p4p . and Ut is the desired target single-qubit opera- 
tion, Eq. UJ). 

We have found numerical solutions involving N = 1, 3, 
and 5 steps for the H, T, and S gates respectively. Ex- 
plicit sequences are shown in Table HP^ where the time 
pulse duration is expressed in terms of the dimensionless 
parameter r = t/to with t$ = (2-7T / 'u; x )[8ir 2 ?iu x /Ry] 1 / 2 . 
Here, Ry = m*e 4 /(2e 2 ?i 2 ) is the effective Rydberg en- 
ergy, m* is the effective mass, and e the dielectric con- 
stant. For GaAs, Ry « 5.93 meV and to « 2.5 ps for 



4 



TABLE I: Pulse sequences for one-bit logic gates. The di- 
mensionless parameters r, z, and r can be tuned through the 



time t and any two of u> y 



and B z 



Hadamard gate 
2ttt z r 



gate 



2-KT 



Phase gate 
2ttt z r 



1.470 0.376 0.158 1.859 4.828 0.022 2.092 0.249 0.121 

3.674 0.102 0.936 1.512 2.803 0.996 

2.443 1.093 0.051 2.123 2.586 0.012 

2.280 0.124 0.916 

1.992 0.224 0.139 




FIG. 3: (Color online) Schematic of our minimal coupling 
model. As the interdot barrier is lowered, the orbitals highest 
in energy and with the greatest overlap will be the first to 
couple. 



u) x = 2.5 meV. Note that the sequences shown in Table U 
are not optimized for experimental efficiency, but merely 
demonstrate that solutions for a universal set do indeed 
exist. We have found many more solutions (none shorter) 
for each of the one-bit gates, including solutions at fixed 
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Z. 



B. Coupled Dots 

To consider two-qubit gates, we now include the inter- 
dot coupling term H coup \. We consider a minimal model 
valid in the limit of weak coupling. Of the three orbitals 
we are considering, the \nm) — 1 02) orbital is both high- 
est in energy and closest to the edge of the dot. Thus, 
as the inter-dot barrier is lowered, the respective 1 02) 
orbitals in each dot will be the first to couple. This is 
schematically depicted in Fig. [3J Our minimal model 
considers the coupling only between these two orbitals. 



1.2 



This leads to a Heisenberg form,-- 

^coupl — J ■ S r , 



(15) 



where S = J2 SS < cl 2s cr ss > c 02 s 
1 02) orbital, and the indices i 
dot respectively. 



is the spin operator of the 
r denote the left and right 



C. Double Qubit Gate 

A two-qubit system is formed from the direct prod- 
uct \Q) ® \Q'), of the states in Eq. ([7]), forming a four- 



TABLE II: Cnot implementation with an always-on intradot 
exchange interaction. Subscripts denote individual qubits, J 
the interdot exchange, and r the pulse duration. 
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Z'2 


7'2 


1.227 


2.133 


0.846 


0.630 


3.280 


0.398 


3.821 


0.615 


1.860 


0.067 


0.663 


0.308 


2.766 


4.094 


0.418 


0.767 


3.897 


0.340 


1.167 


3.540 


0.017 


0.298 


0.852 


0.952 


1.591 


3.242 


1.695 


0.370 


2.362 


0.237 


2.148 


3.031 


2.177 


0.559 


2.648 


0.354 


1.560 


1.714 


3.091 


0.077 


4.812 


0.083 


2.255 


1.889 


1.536 


0.222 


2.032 


0.645 


1.981 


3.796 


21.501 


0.453 


11.516 


0.157 



dimensional computational space. The states \Q) <£> \Q') 
are, in fact, S = 1, S z = —1 eigenstates. Unfortu- 
nately, the spin subspace and the computational sub- 
space are not identical; the six spin (three for each dot) 
S = 1, S z — — 1 subspace is nine-dimensional, four of 
which constitute our \Q) <g> \Q') computational space. 
Thus, our implementation of cnot, as that in Ref. [T^, 
involves transient excursions outside the computational 
space; nevertheless, our sequences are designed such that 
the final gate operation is unitary and returns to the four- 
dimensional computational space. We require the final 
state to be such that the cnot truth table be satisfied, 
up to single-qubit operations^ An explicit implementa- 
tion is given in Table [TT1 Six parameters are required to 
describe a pulse: the pulse duration r, two (r, z) param- 
eters per dot to describe each qubit, and a dimensionless 
exchange coupling J — JtTi/r describing the coupling. 
As shown in Table UH a nine step solution is the small- 
est we have been able to find^S (If it were possible to 
turn off intradot exchange, then a three-pulse cnot is 
achievable^) 



V. DECOHERENCE AND DISSIPATION 

Environmental influences can be of two distinct types: 
Slow variations in the electromagnetic environment 
merely lead to adiabatic changes in the pseudofield b 
and thus to unitary errors that typically average out 
over the length of a pulsed We look first to fast, non- 
adiabatic environmental influences that can lead to non- 
unitary errors — i.e., decoherence — followed by a discus- 
sion on adiabatic influences which lead to gate errors. 



A. Nonadiabatic Influences 

Assuming the environment does not change the num- 
ber of particles on the dot, and that the bath couples 
only to single particles in the dot, then a general model 
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of system-bath coupling is given by 

^SB = ^ -B n ms C^, m , s ,C nms , (16) 

where the sum is over all repeated indices. B™™ s s is 
a set of arbitrary operators which describe the reservoir 
and all relevant coupling constants. 

At time t, the full state- vector of the system |$(t)) = 
J2q \Q) ® IXq(*)) (Q = 0' 1) 1S a non-separable state, 
where the states \xc}(t)) are reservoir states including all 
time-dependent coefficients. Matrix elements of Eq. (|T6"|) . 

^S Q B Q = £<Q'l4'mV<W.|Q) A™'( X ',X), (17) 

where A™Z s '(x',x) = (x'\KZ s '\x) are straightfor- 
wardly calculated.— 

Within the Born-Markov approximation, and using the 
definitions Eq. 0, the relaxation T\ and dephasing T v 
times are given by^ 2 . (1/T 2 = 1/(2T X ) + 1 /T v ) 

Y x ~\ H ^\ 2 = h 15/10 ~ 5/11,2 ' (18a) 

i- ~ |# s °° - i^| 2 = ± |2M a - [Shi + Sh )\ 2 , (18b) 

where 5h m — ^4q™{ — A^mi To the extent that the 
Born-Markov approximation is valid, 23 Eq. (| 18[) states 
that relaxation and dephasing within the computational 
space are negligible to leading order for all environmen- 
tal couplings which are either purely charge or purely 
spin in character. The former has B™™ a s — (5 ss /£?™ n ™ 
in Eq. (fTrJ)) and consequently Sh m = 0, whereas the lat- 
ter has B™™ 3 s = S m ' m S n ' n Bg and consequently Sho — 
Shi = 5h,2 7^ 0. For both these cases, dephasing and re- 
laxation vanish within the Born-Markov approximation. 
(Neither of these is applicable for hyperfine environments 
which depend on both spin and charge.) 



for example, for the values r w 0.8 and z — > 0. For 
these values of the confinement, we find \db/dz\ — > 
and \db/dr\ ~ e 2 /(16\/27rA), where A 2 = K/(mu) x ). For 
GaAs material parameters, with lj x = \ meV, this gives 
\db/dr\ « 85 /ieV — at least an order of magnitude smaller 
than typical pscudofield magnitudes. With these, we can 
estimate corrections to the adiabatic limit^l expressed 
as a leakage time given by 7i ca k ~ E 4 /(h 2 r 3 \db/dr\ 2 ), 
where E ~ 1 meV is the excitation energy to states out- 
side the qubit space, and f = dr/dt ~ 10 GHz is the rate 
of typical gate operation. For these parameter values, we 
find Ti ca k ~ 200 /j,s, and a leakage probability of only 

VI. INITIALIZATION AND MEASUREMENT 

With regard to state initialization, we note that if the 
qubit states are the lowest-energy states, then prepara- 
tion becomes merely a matter of thermalization. In this 
case, potentials other than elliptic may well prove useful. 9 
Ideally, a system where the two qubit states are the two 
lowest-energy states would be beneficial not only for state 
preparation, but also for a more general deterrent to dis- 
sipation, especially relaxation to states outside the com- 
putational basis. 

Finally, with regard to measurement, we note that the 
two states in Eq. (J7J) are not degenerate. Thus, a de- 
structive measurement is possible by detecting whether 
a fourth electron resonantly tunnels onto the (three- 
particle) dot; similarly to the single-shot readout of in- 
dividual quantum dot spin^ the gates may be pulsed 
such that an additional electron can tunnel onto the dot 
only if it is in the higher-energy qubit state. In fact, 
since the tunnel barriers as well as the confinement it- 
self is determined (and controlled) by the applied elec- 
trostatic potential, the universal set of gates described 
above as well as detection may be accomplished using 
already-existing^ 2 ^ experimental techniques. 



B. Adiabatic Influences 



Regarding adiabatic (unitary) influences, which do not 
cause decoherence, these can be minimized by choosing 
settings for the confinement potential such that db/dr 
and db/dz are small in magnitude. This is the case, 
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